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ABSTRACT 


This report demonstrates the feasibility of obtaining measure- 
ments of Lagrangian turbulence at stratospheric altitudes by using 
the METRAC System to track constant -level balloons. The basis for 
current estimates of diffusion coefficients are reviewed and it is 
pointed out that insufficient data is available upon which to base 
reliable estimates of vertical diffusion coefficients. It is 
concluded that diffusion coefficients could be directly obtained 
from Lagrangian turbulence measurements. The METRAC balloon tracking 
system is shown to possess the necessary precision in order to resolve 
the response of constant- level balloons to turbulence at stratospheric 
altitudes. A small sample of data recorded from a tropospheric 
tetroon flight tracked by the METRAC System is analyzed to obtain 
estimates of small-scale three-dimensional diffusion coefficients. 

It is recommended that this technique be employed to establish a 
climatology of diffusion coefficients and to ascertain the variation 
of these coefficients with altitude, season and latitude. 
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1. INTRODUCTION 

The diffusion of pollutants introduced into the stratosphere by 
supersonic transports has been the subject of an extensive research eflort 
coordinated and sponsored by the Department of Transportation. This 
research program known as CIAP (Climatic Impact Assessment Program) has 
brought the resources of several government agencies, universities, and 
private institutions to bear on an evaluation of the danger of a depletion 
of the earth's ozone shield as a result of the operation of a fleet of 
supersonic transport aircraft flying at altitudes from 18 to 21 km in the 
stratosphere. One of the central objectives in evaluating this danger has 
been the development of realistic models which incorporate both the effects 
of turbulent diffusion and the effects of chemical transformations. These 
models can be one -dimensional, two-dim nsional or three-dimensional. The 
two-dimensional and three-dimensional models are, of course, more realistic 
than one -dimensional models but they require more time for development and 
consume much more time on the computer. In all of these models diffusion 
coefficients are required to account for the turbulent spread of effluents. 
In order to make the problem tractable it is necessary to consider several 
different time and space scales. The smallest scale of interest is the 
scale of the wake of the stratospheric transport. The wake itself can be 
subdivided into jet, vortex and dispersion regimes (Poppoff, Far low and 
Anderson, 1974). In the dispersion regime a transition occurs from aircraft 
induced dispersion to dispersion dominated by the turbulence of the natural 
stratosphere. This occurs between 100 and 1000 sec after passage of the 
aircraft. At ten minutes after passage of the SST the width of the wake is 
about 500 meters and its depth is about 16 meters (Hoshizaki, Anderson and 
Conti, 1972). 

Hie smallest atmospheric scales of interest in the dispersion of the 
aircraft wake are on the order of five minutes and 300 meters. Beyond 
these scales, diffusion is controlled by stratospheric wind fields and the 
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wake spreads and mixes with other wakes. After several weeks the problem 
becomes global in scope. 

This report is directed toward an examination of a new approach to 
the in situ measurement of small-scale diffusion coefficients in the 
stratosphere. This new approach is to accurately track constant -level 
balloons in three dimensions using the METRAC System and to deduce diffu- 
sion coefficients from an analysis of Lagrangian turbulence statistics. 

The METRAC System is a radio-location system which employs the Doppler 
principle in order to accurately position an expendable radio transmitter. 
Direct measurement of small-scale vertical motions in the stratosphere 
using METRAC may provide the only means for direct determination of vertical 
diffusion coefficients. 

The nature of stratospheric turbulence and the basis for current 
estimates of stratospheric diffusion coefficients are presented in Chapter 
2 and Chapter 3. A review of balloon techniques for obtaining turbulence 
measurements and diffusion coefficients is presented in Chapter 4, Chapter 

5 contains a description of the METRAC balloon-tracking system and Chapter 

6 contains an analysis of a sample of data from a tetroon tracked in the 
troposphere by the METRAC system. Some considerations of the application 
of the METRAC system for tracking constant- level balloons in the strato- 
sphere are discussed in Chapter 7. Conclusions and recommendations are 
given in Chapter 8. 

2. NATURE OF STRATOSPHERIC TURBULENCE 

In order to understand the scale-dependent spread of atmospheric 
pollutants in the stratosphere it is necessary to review the nature of 
turbulence in the stratosphere. Classical turbulence theory has been 
developed for statistically homogeneous, isotropic fields of turbulence. 
Although turbulence approaches the homogeneous isotropic state in the lower 
atmosphere, it can be isotropic only on the smallest scales in the strato- 
sphere (A ~ 10 m, t - 100 sec). This is primarily due to the stable thermal 
stratification which inhibits vertical motion. Thus it is common experience 
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that stratospheric effluents diffuse more vapidly in the horizontal than 
they do in the vertical. The nature of turbulence in a stably stratified 
fluid under the influence of the gravitational force is intermittent in 
time and inhomogeneous in space. Dynamically, turbulence develops in local 
zones of strong shear. The Richardson number is commonly used as a measure 
of the dynamic stability of a stratified shear flow. When the local 
gradient Richardson number 


(0 is potential temperature) 

is reduced below 1/4, an instability and transition to turbulence may occur. 
The resulting turbulence tends to destroy gradients within the turbulent 
layer and concentrate them at the boundary. Thus, statically stable zones 
are often found to contain much small-scale layered structure; with weakly 
turbulent zones separated by stable laminae. Gradients are concentrated in 
the stable laminae which occasionally break down and become turbulent. 

This structure is typically found in the stably stratified ocean on a very 
small scale. In the thermocline region stable laminae may be on the order 
of a few tens of cm thick [Woods (1969)]. In the stratosphere a layered 
structure also exists but it has not been possible to observe this 
structure in the same detail as in the ocean. 

In order to model the spread of stratospheric pollutants in three 
dimensions it is necessary to parameterize the effect of turbulent diffusion. 
Because of the anisotropic nature of stratospheric turbulence it is important 
to consider separately vertical and horizontal diffusion. If it were 
possible to record an infinitely long time series of stratospheric winds, 
analysis would probably yield a spectral peak of the horizontal wind spectrum 






in the large-scale (synoptic and global scale) end of the spectrum with a 
decrease in energy toward smaller scales and a second peak near the micro- 
scale* This is shown schematically in Fig* 1 adapted from Vinnichenko 
(1970)* The existence of the meso-scale gap observed by Van Der Hoven 
(1957) in the boundary layer has been explained on theoretical grounds by 
Fiedler and Panofsky (1970). Observations are not available to confirm 
its existence in the free atmosphere. If direct vertical velocity measure- 
ments were available in the free atmosphere in a time series, an analysis 
would reveal a much different spectral distribution than for the horizontal 
components of velocity. In the asymptotic limit of the smallest scales the 
horizontal and vertical velocity spectra would appear quite similar. At 
larger scales, however, the vertical velocity spectrum would contain a 
negligible amount of energy compared to the horizontal component shown 
schematically in Fig* 1, In fact, the peak of vertical velocity spectra 
would occur on a time scale of only a few minutes [Panofsky (1969)]. 

The different distributions of vertical and horizontal turbulent 
kinetic energy are responsible for the different rates of spread of pollu- 
tants in vertical and horizontal directions. It is generally understood 
that maximum spread or diffusion occurs on or above the scale of maximum 
turbulent energy. This view J s consistent with scale dependent diffusion 
coefficients which increase with scale through regions of strong turbulent 
energy and which approach constant values at scales beyond which there is 
little additional turbulent kinetic energy. Thus, schematically as shown 

in Fig. 2, adapted from Bauer (1974), the lateral diffusion coefficient 

- - - - - ~ 2 2 - 1 -- — 
increases with scale from 10 cm s at a 10 second time scale to 

10^ cm^ above 10^ sec. The vertical diffusion coefficient increases 

2 2 -1 4 2 -1 3 3 

from 10 cm s at 10 sec to only 10 cm s above 10 sec. Above 10 

sec the vertical diffusion coefficient approaches its asymptotic upper 

limit. 
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Figure 1, Schematic representation of frequency 
spectra of horizontal and vertical 
velocities in tne free atmosphere 
(Zonal wind spectrum adopted from 
Vinnichenko, 1970). 



Dispersion Distance (km) 



Figure 2. Schematic representation of the variation 

of lateral and vertical diffusion coefficients 
j with diffusion time (after Bauer, 1974). 
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3. REVIEW OF STRATOSPHERIC TURBULENCE MEASUREMENTS AND THE BASIS FOR 

ESTIMATES OF STRATOSPHERIC DIFFUSION COEFFICIENTS 

In the course of the CLAP program an effort has been made to review 
the sources of stratospheric turbulence measurements and to infer strato- 
spheric diffusion coefficien s. The results of such a study have been 
reported by Bauer (1974). On the small scale, turbulence measurements have 
been made and diffusion coefficients inferred from aircraft and from the 
dispersal of smoke puffs. On the global scale, diffusion coefficients 
have been inferred from measurements of radioactive material injected into 
the stratosphere from nuclear explosions and from global heat flux data. 
Below are reviewed the principal techniques for obtaining estimates of 
diffusion coefficients. Vertical diffusion coefficients are considered in 
more detail since their basis is much less certain than that of lateral 
diffusion coefficients. 

3.1 AIRCRAFT MEASUREMENTS 

Between 1964 and 1968 the Air Force recorded approximately 48 
hours of turbulence data from specially instrumented U-2 aircraft flown 
between 15 and 21 km. These data were taken for the Air Force HICAT prograiu 
and are the subject of several technical reports [Crooks, et al. (1967, 
1968) and Ashburn, et al. (1968, 1969, 1970)] written for the Wright 
Patterson Air Force Flight Dynamics Laboratory. 

Instrumentation aboard the HICAT U-2 aircraft provided measurement 
of all three components of velocity as well as temperature. The accuracy 
of these measurements is discussed by Crooks, et al. (1968) but a complete 
error analysis has never been made. 


I' 


Lilly, Waco and Adelfang (1974) have recently reviewed the HICAT 
data in order to deduce estimates of stratospheric diffusion coefficients. 
The method used is as follows. First, the dissipation (c) is estimated 
from composite spectra assuming an inertial subrange form for the spectral 
density function. An assumption that the flux Richardson number should 
be near 1/4 leads to the conclusion that buoyancy or heat flux is equal to 
one third of the dissipation. Then, if the usual Fickian relationship is 
assumed to hold between heat flux find temperature gradient, it follows that 
the diffusion coefficient for heat is 


where 




e az . 


(3.2) 


N is the Brunt-Vaisala frequency, and 0 is potential temperature. Now an 

o -4 ~2 

isothermal stratosphere is assumed [with T = 210K, N = 4.6 x 10 sec ] 

and the heat flux can be evaluated from the dissipation. This leads to an 

^,:direct estimate of diffusion coefficients from aircraft measurements of the 
2 

\ locity variance v ( k) . The values of the estimated diffusion coefficients 

3* 2 -1 4 2-1 

range from 4.8 x 10 cm s over water to 1.2 x 10 cm s over high mountains 

Since these estimates are only applicable to regions of intermittent 

turbulence^ a more representative diffusion coefficient was estimated by 

multiplying the above numbers by their probability of occurrence. Hie result 

is an estimate of an effective diffusion coefficient which ranges from 
2-1 2-1 

100 cm s over land to 640 cm s over high mountains. These estimates 
were made for turbulence on the scale of 610 meters which corresponds to 
a few secondS|On the diffusion time scale of Bauer's diagram. 
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Some addidional turbulence data is available from stratospheric 
flights of Project Coldscan reported by MacPherson and Morrissey (1969, 

1970) and from the Colorado LEE WAVE Experiment reported by Lilly, et al. 
(1971), and Lilly and Lester (1974). These data are not nearly as complete 
as the HICAT data reported above. 

3.2 DISPERSION OF SMOKE PUFFS 

Another source of small-scale turbulence data is provided from 
observations of the dispersion of smoke puf. c s and trails. This work is 
most significant since small-scale diffusion may be more naturally assessed 
in a Lagrangian context. The dispersion of smoke puffs was observed as 
part of the "High Altitude Dust Diffusion Project” at Holloman Air Force 
Base. Results of twenty experiments at altitudes ranging from 7.2 to 19.3 
kilometers were reported by Kellogg (1956). Smoke puffs were created from 
the explosion of vials of titanium tetrachloride and water. These vials 
were attached to a train of rising balloons and were exploded by the closing 
of a baroswitch contact at pre-selected altitudes. 

After the initial disturbance created by the explosion died 
down, the resultant smoke puff acted as a tracer and was tracked by a 
network of three photo theodolites . Later the growth of the smoke puff was 
deduced from a detailed hand analysis of photographic data. Specif ically , 
the cloud edges were observed and the diameter of the cloud was determined 
as a function of time. After about six minutes the cloud became so tenuous 
that it was difficult to reliably measure the cloud diameter. In order to 
separate the effect of diffusion from the effect of wind shear the data 
was analyzed to determine the minimum diameter at each time. The average 
values of the minimum diameter are presented in Fig. 4 of Kellogg (1956). 

The puffs typically grew from 15 meters to 120 meters in six minutes. This 
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growth is consistent with a diffusion coefficient of 5 x 10 cm s for 
diffusion times of five minutes. There is no attempt in the analysis to 
differentiate between horizontal and vertical diffusion coefficients. 

3.3 RADIOACTIVE TRACERS 

There are three categories of radioactive material from nuclear 
explosions. These categories are radioactive products of the fission 
process, radioactive isotopes activated by neutron bombardment and unspent 
radioactive bomb material. A common tracer in the first category is 
strontium-90 which has been the subject of considerable study because of 
its long half-life and its potential health hazard. Carbon-14 is produced 
by neutron activation. 

The commonly measured radioactive isotopes include strontium-90, 
carbon-14, tungsten-185, plutonium-238, cadmium-109, and rhodium-102. 

The first two isotopes are produced by every nuclear test v;hereas tungsten- 
185, cadmium-109, and rhodium-102 are special tracers which can be identi- 
fied with certain tests. Plutonium-238 was injected into the atmosphere 
as a result of the burn-up of a navigational satellite with a radioactive 
power source. 

Tracer data has been the subject of extensive analysis by a number 
of authors. The main features of the results are summarized by List and 
Telegadas (1969). Two periods have received the most attention. These 
are the periods of test moratoriums which were 1959-1961 and after January 
1963. 


Results of the analysis of tracer data are incomplete. However, 
some obvious conclusions can be drawn from available data. It is clear that 
the spread of radioactive material is dependent upon altitude, latitude and 
season. Residence times give an indication of the rate of vertical trans- 
port of radioactive tracers. In the lower equatorial stratosphere tungsten-185 
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tracers produced in 1958 exhibited a residence time of eight months. Tracer 
data, however, did not become available routinely above 20 km until 
October 1964, so that a reliable estimate of residence time for material 
injected into the upper stratosphere is not available. Nevertheless, 
Peterson (1970) has given an estimate of two years for the residence time 
there. Russian tests in the polar regions in 1958 produced tracers which 
were confined below 15 km from which it was concluded that residence time 
for radioactive materials injected into the lower polar stratosphere is 
about five months (Peterson, 1970). 

One of the difficulties in interpreting the results of tracer 
dispersion is that the tracer spreads due to the combined influence of . 
large-scale circulations and diffusive mixing. List and Telegadas (1969) 
tentatively conclude from their analysis of tracer data that there are 
regions of the stratosphere which are dominated by advective transport 
and other regions dominated by mixing. The altitude range from 18 to 23 
km, and the latitude belt from 25N to 2.5S is thought to be dominated by 
mixing. Hie entire summer stratosphere below 40 km is similarly thought 
to be dominated by mixing processes. List and Telegadas (1969) also con- 
clude that at least between 25° and 70° latitude there is a strong mean 
descending motion in the winter stratosphere which dominates vertical 
transport. 

There has been some concern expressed recently by Johnston, et al. 
(1975), that residence times deduced from particulate radioactive tracers 
are seriously underestimated. These authors advocate the use of carbon-14 
as a tracer more representative of atmospheric motions. Carbon-14, of 
course, is found in CO^ which is a gaseous tracer. The possibility that 
radioactive particulates will settle under the influence of gravity has been 
explored by Telegadas and List (1969). . Settling rates are dependent upon 
the radius of the particle and the atmospheric density at the altitude 
being considered. Telegadas and List conclude that the settling rates are 
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too small to be important for the size of particles of interest below 25 to 

14 

30 km. Their conclusion is supported by an analysis of ratios of C to 
90 

Sr .- For example, between 1959 and 1962 the ratio of these two tracers is 

about constant with time at selected altitudes and latitudes. The fact 

that strontium-90 decreases with time more rapidly than carbon- 14 during 

the second test moratorium is attributed to a difference in the initial 

vertical distributions of the two tracers. In this connection it is believed 

14 

that Russian tests in 1961 produced much higher concentrations of C relative 
90 

to Sr at altitudes above 25 km. 

In view of the analysis of Telegadas and List (1969) and the radio- 
active tracer data reviewed above it is concluded here that the season, and 
location of stratospheric injection are the determining factors in explaining 
the differences in residence time between C^ r and particulate radioactive 
tracers pointed out by Johnston, et al. (1975). . 


3.4 ESTIMATES BASED ON HEAT FLUX DATA 

Another source for diffusion coefficients is the work of Reed and 
German (1965). These authors compute seasonal global diffusion coefficients 
based on heat flux data and employing mixing length assumptions. Reed and 
German use the relation between vertical diffusion coefficient and 
lateral diffusion coefficient Kyy given by 

Hz " (<* 2 -+ o ,2 )^pi - . _ . , . . . (3.3) 

_ 2 

where a is the zonal mean slope of the mixing surface and a* is its 

variance. The lateral diffusion coefficient Kyy is determined as a function 

of latitude from the variance of the meridional velocity component. From 

the heat flux and temperature gradient data it is possible to deduce the 
- 

variation of a with latitude but a* is unknown. Reed and German use a 

3 2-1 

representative equatorial JC value of 10 cm s deduced from the vertical 




spread of radioactive tracers to ob tain a value of cr at the equator where 
ct = 0. Then, "Not knowing how o' varies with latitude and height, we 
assume that it remains constant/ 1 In this manner and from the previously 
deduced variation of a with latitude a latitudinal variation of is 
generated. 

There are several reasons to doubt the validity of the large- 
scale thermal diffusion coefficients One reason is that the computation 

of is not direct but rather is adjusted to agree with a. vertical diffu- 
sivity obtained from tracer studies. However, even if the approach of Reed 
and German gives valid vertical thermal diffusion coefficients it is not 
clear what relationship these thermal diffusion coefficients will have to 
the desired diffusion coefficients of passive additives. Diffusion coeffi- 
cients are defined by an assumed linear relationship between flux and gradient. 
On the global scale, mixing surfaces and surfaces of constant potential 
temperature are inclined to the horizontal and do not share the same slope. 
Because of this situation, the gradients and fluxes may be quite different 
for heat than for a passive contaminant. The common assumption that diffu- 
sion coefficients should be nearly the same for different transported 
quantities is based on the hypothesis that the physical mechanism causing 
transport is common to the quantities being transported. That this is 
approximately true for small-scale turbulent transport is well documented. 

For large-scale transport it is not true that the physical mechanism respon- 
sible for vertical heat flux is common with mass flux. 

3.5 DISCUSSION 

The spread of material injected into the stratosphere is affected 

by the entire spectrum of atmospheric motions. For the horizontal spread 

the largest scale motions dominate the transport process and global diffusion 

9 2-1 

rates have been estimated to be on the order of 10 cm s . There are 
several ways in which this estimate can be supported. What heeds clarifica- 
tion is how this global diffusivity is approached from the small-scale end 







with increasing diffusion time. 

Estimates of vertical diffusivities are in general not well 
supported because of insufficient data. Perhaps the primary reason for 
this state of affairs is the lack of direct measurement of atmospheric 
vertical velocities. We do know that mean vertical velocities are extremely 
small except on the smallest scales. If it is true that the spectrum of 
vertical velocity has its peak on the scale of a few minutes, it seems reason- 
able that this scale will dominate the vertical transport process. Unfortun- 
ately, there is just not enough data available to ascertain which scales of 
atmospheric motion do in fact dominate the vertical transport process. 

The available information for estimating vertical diffusion coeffi- 
cients in the stratosphere can be summarized as follows. Lilly, etal, (1974), 

have analyzed HICAT aircraft turbulence data and concluded that typical 

2-1 

vertical diffusivities are on the order of 100 cm s . The estimates of 

Lilly, et al. (1974) are considered valid for diffusion times on the order 

of 5-10 seconds. Kellogg (1956) has presented an analysis of the spread of 

smoke puffs from 18 stratospheric trials. His data is consistent with a 

diffusion coefficient of roughly lO^cm^ s \ for diffusion times on the order 

of five minutes. Unfortunately it is not possible to separate vertical from 

horizontal transport in Kellogg* s data so that the number quoted above may 

be an overestimate of vertical diffusivity. Finally, the most definitive 

data for basing estimates of diffusion coefficients comes from analyzing the 

spread of radioactive tracers. Vertical diffusion coefficients on the order 
3 2- -l 

of 1-5 x 10 cm*' s~ (Machta and Telega'das, 1972) are consistent with the , _ . 

vertical transport of radioactive tracers. It has not beer, possible, with 
available data to obtain variations of vertical diffusion coefficients with 
latitude, altitude and season. The latitudinal variation of K 72 presented 
by Reed and German (1965) is based. on several questionable assumptions and 
should not be used uncritically. 




4, balloon techniques for measuring turbulence act determining diffusion 

COEFFICIENTS 


Turbulence theory is developed in two frames of reference: Eulerian 

and Lagrangian. In tho Eulerian frame a sensor is fixed in space and 
the time history of the fluctuating quantity being measured at that point 
is recorded* Time series analysis can be employed to analyze the frequency 
content of the turbulence and the mean velocity can be used according to 
Taylor's hypothesis to relate temporal to spatial variations. For obvious 
reasons the Eulerian frame is the preferred frame of reference for all 
surface boundary layer studies and extensive Eulerian turbulence data is 
available from tower mounted instruments . The Lagrangian frame follows 
the motion of elements or parcels of fluid. Turbulence measured in this 
moving frame is referred to as Lagrangian turbulence. Lagrangian turbu- 
lence is more basic to diffusion theory since it is in the Lagrangian 
frame that material actually diffuses. 

Since most atmospheric measurements near the surface are taken in the 
Eulerian frame, considerable effort has been made by micrometeorologists to 
relate Eulerian turbulence statistics to Lagrangian turbulence statistics. 

Above the surface boundary layer it is much more difficult to collect meaningful 
in situ data. The lack of fixed platforms in the free atmosphere dictates 
the search for alternative techniques for measuring turbulence. The 
analysis of the spread of smoke puffs and radioactive tracers has already 
been discussed. These tracers offer the advantage of Lagrangian measure- 
ments but they do not readily lend themselves to quantitative analysis. 

Two techniques which do lend themselves to quantitative analysis are measure- 
ments from aircraft and rising balloons. Neither of these techniques 
yields an Eulerian or Lagrangian view of turbulence. 

A far more powerful balloon technique for measuring turbulence in 
the free atmosphere has been developed in the past fifteen years. This 
method is to track r.onstant" level balloons. Data so obtained can be readily 
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interpreted in terms of Lagrangian turbulence theory and the results are 
immediately applicable to diffusion theory. 

Rising or floating balloons can be used to study air motion on a 
variety of scales. When attempting to extract turbulence information 
from balloon motion, two considerations are of vital importance. The 
accuracy with which the balloon can be positioned in space is the subject 
of Section 4.1. The degree to which the balloon motion provides an accurate 
indication of atmospheric motion is discussed in Section 4.2. The 
remainder of the chapter is concerned with two techniques for deducing 
diffusion coefficients from balloon measurements of turbulence. 

4.1 REVIEW OF BALLOON TRACKING SYSTEM 

All balloon tracking systems currently available use one or a combina- 
tion of the following techniques to obtain position or wind information: 

1) Azimuth and elevation angles 

2) Passive ranging (reflection) 

3) Active ranging (transponders) 

4) Thermodynamic height evaluation 

5) Navigational aids 

6) Doppler frequency measurements 

The following paragraphs will sutrenarize the wind finding or position- 
ing capabilities of eight types of balloon tracking systems. 

o Single Theodolite 

Single theodolite measurements have been shown to yield reasonable 
wind profiles in the lowest few kilometers. Angular accuracy with optical 
theodolites can approach 0.01°, but observations are limited by cloud cover. 
Radio-theodolites have stated accuracies of only about 0.1° but are not 





limited by cloud cover or other visual obstructions. Accuracy for both 
types of systems is angle dependent. The necessary assumption of a constant 
balloon ascent rate implies large wind errors for both very high and low 
elevation angles. Deviations in the balloon's ascent rate are translated 
into horizontal wind speed errors. In summary, the single theodolite system 
is of practical use in obtaining limited resolution vertical profiles of 
horizontal wind in the lowest few kilometers of the atmosphere. 



Multiple theodolite tracking of balloons has been used to obtain more 
detailed wind profiles than are possible with single theodolite observations 
(e.g., Ackerman, 1974). The balloon position, given geometrically by the inter- 
section of lines from the theodolites, is independent of any assumption of 
the balloon's vertical motion or any form of atmospheric equilibrium. Since 
there are always some errors in the observations, the most common method for 
solution is to follow Thyer (1962) who positions the solution along the line 
connecting the points of closest approach between the rays defined by the 
theodolite angles. Radio-theodolites allow observations during cloudy days 
and, with large baseline separations, up to stratospheric heights. Nelson 
(1973) has analyzed data collected from GMD-2 and WBRT-60's operated in. a 
triangular array approximately 50 km on a side. Even though equipment 
accuracy specifications implied an uncertainty in positioning of 40 to 80 m 
for any two rays of total length 50 to 100 km, observed errors were an 
order of magnitude larger. Nelson attributed these errors to a combination 
of antenna alignment problems and to overly optimistic equipment accuracy 
specifications. The former will yield a systematic or absolute error while 
the latter will produce a random positioning error. The random error is of 
greater significance if one wants to look at turbulence structure. Smaller 
baselines between radio-theodolites will reduce the effects of errors caused 
by both systematic and random error components for times near launch. However, 
errors in balloon positioning at large distances and upper tropospheric or 
stratospheric heights will increase. 
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• Angle Positioning with Thermodynamic Heights 

This method has been used historically by the Weather Service and others 
to determine air motion with the radiosonde (rawinsonde) . Knowledge of the 
thermodynamic state parameters allows integration of the hypsometric equation 
to determine height which together with elevation and azimuth angles yields 
the balloon position. Errors in wind will arise from errors in the measured 
angles as well as errors in height due to inaccuracies in the numerical 
integration of the hypsometric equation or departures of the atmosphere 
from hydrostatic balance. The GMD series of tracking systems have been 
the standard in radiosonde tracking for about a quarter century. Commonly 
accepted rms positioning errors are 0.05° in azimuth and elevation angles 
for elevation angles above 15°, though Danielsen and Duquet (1967) show 
evidence of occasionally larger and nonrandom angular deviations from the 
true position which can persist for several minutes. The assumption of a 
uniform ascent rate between standard pressure contact levels can also lead 
to errors in determining statistics of the horizontal wind field. All in 
all, errors in the horizontal wind of +1 m are typical over one minute 
Intervals when angular data is sampled ten times each minute. Tracking 
angles below 15° produce larger errors due to ground reflection and refrac- 
tive index variations. Clearly, this technique is unsuitable for tracking 
const ant -level balloons. 

• Angle Positioning with Passive Ranging 

Passive ranging of a balloon sensor has been demonstrated to be one 
of the more accurate methods of determining small-scale atmospheric motion. 
Sophisticated radars like the FPS-16 have been able to resolve motion on 
scales smaller than the natural induced oscillations of rising standard 
meteorological balloons. Equipment specifications give rms accuracies of 
.01° in azimuth and elevation and five meters in slant range. However, 

Scoggins and Armendariz (1969) discuss some of the practical limitations in 
accuracy of the FPS-16 radar system. Operator adjustments and data smoothing 
procedures are very important for optimizing data quality. DeMandel and 
Krivo (1972) have determined the frequency at which most of the variance 
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in the observations can be attributed to atmospheric motion rather than 
radar noise. They found that structure with a vertical scale of 10 
meters was resolvable at heights of 5 km, however, at 15 km, only structure 
with a vertical scale over 100 meters is observable. Also, as is the case 
with all angle positioning devices, errors are substantially increased at 
low elevation angles. This technique has been used to track constant -level 
balloons as well as rising balloons. 

o Angle Positioning with Active Ranging 

Another method of balloon positioning is to attach an active device 
(transponder) to the balloon for ranging purposes. This technique used 
in conjunction with an M-33 radar has allowed accurate tracking of low level 
tetroons to distances of 105 km by research teams at the National Reactor 
Testing Station (NRTS) in Idaho Falls, Idaho (e.g., Angell, et al., 1968). 
Estimates of errors associated with this system have been made by simultaneously 
tracking a single tetroon with two M-33 radar systems. These data give results 
of rms system errors of about 30 meters in range and nearly 0.10° in azimuth 
and elevation angles. These numbers translate into instrument errors of 
about 2.5 m s" 1 in the horizontal wind and 1ms" 1 in the vertical air motion 
over periods of 30 seconds. Much of the M-33/tetroon data is smoothed over 
time periods of two to three minutes. 

• Navigational Aids Positioning 

Within the last ten years a new type of wind measuring system has 
been developed which makes use of already existing VLF and LF frequencies 
available for general navigation. Principal Navigational Aids (Navaids) 
used so far include OMEGA (10-14 KHz) and Lor an (around 100 KHz) trans- 
missions. Many papers [e.g., Acheson (1970), Beukers (1972), Beukers (1973)] 
have been published describing Navaid systems and associated accuracies. 

The Omega system can be used over extensive areas of the globe but is 
limited to accuracies of horizontal wind of around 1ms 1 over four minute 
periods and 3ms 1 for one minute averages (Lally, 1972). Loran systems 
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have accuracies of around 0,5 m for one minute intervals (Beukers, 1975) 
but geographical coverage is limited. Vertical motion and height are 
determined thermodynamically as with a radiosonde. The greatest advantage 
of the Navaid system is for tracking balloons at sea, as it does not depend 
upon the expensive stable platforms required by r.ll other systems. 


• Doppler Radar 

Doppler radar systems can be used to determine fine-scale velocity 
structure indicated by balloon sensors, A feature of this type of system 
is that the radial velocity of the target is measured directly. Doppler 
systems such as the one discussed by Lhermitte (1967) have the capability 
of measuring the instantaneous radial velocity of spherical balloons to 
about 0.5 m s A problem with using a Doppler radar to look at small-scale 
velocity structure is that a single system can only detect one component 
of the velocity field. Three systems separated in space and looking at 
the same target are necessary to obtain the three-dimensional velocity 
field. 

• Doppler Positioning 

Recently, a new type of Doppler positioning system which has the 
capability of looking at very fine scale atmospheric velocity structure 
has been demonstrated. This system, the METRAC™ positioning system, will 
be described in more detail in the next chapter. Here it will suffice to 
say that the METRAC system eliminates most of the problems inherent in 
angle measui ^ systems. Problems of refractive ray bending are minimized, 
and areal separation of system components allow very accurate positioning 
over extended distances although, as will be shown later, this accuracy 
is dependent upon system geometry. 

4.2 BALLOON AERODYNAMICS 

An evaluation of the accuracy of fine-scale measurements of atmospheric 
winds from balloon systems requires consideration of the aerodynamics of 
balloons . 
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• Self-Induced Oscillation of Rising Balloons 

It has been known for decades that rising rubber balloons undergo self- 
induced motions even in still air. Although some investigators have attri- 
buted these motions to deformation of the balloons, it is now known that 
rigid plastic balloons also suffer self-induced oscillations. Murrow and 
Henry (1964) have performed a number of experiments on several prototype 
balloons which elucidate the extent of the problem. Their experiments 
consisted of tracking ascending balloons in Hangar Number One at the Lakehurst, 
New Jersey, Naval Air Station. The balloons tested included standard 
radiosonde balloons, 2-meter mylar spherical balloons, streamlined balloons 
and roughened balloons. All balloons tested exhibited self-induced oscilla- 
tions although some reduction in amplitude of oscillation was achieved by 
significantly roughening the surface of balloons. 

The theory behind the self-induced oscillation is far from complete 
but it is generally understood that the oscillations are the result of 
periodic x and y directed aerodynamic lift forces which result from vortex 
shedding in the unstable wake of the balloon. According to Fichtl (1971), 

TO) — 

as long as the Reynolds number (Re = — ; where w is the average rate of 
vertical ascent, D is the diameter of the balloon and v is kinematic viscosity) 
is below a critical value of 2.5 x 10^, a rising spherical balloon will 
exhibit an orderly spiral motion with a vertical wavelength of close to 
twelve balloon diameters. However, for super-critical Reynolds numbers, 
the balloon^ self-induced motion becomes erratic. For a spherical balloon 
of two meter diameter with an average vertical velocity of five meters per 
second in the standard atmosphere, the Reynolds number will be super- critical 
below 11 km and subcritical above. 

The Jimsphere is a spherical super-pressure balloon (Scoggins, 1965), 
the surface of which has been covered by roughness elements. The Jimsphere 
has been developed by NASA and used in conjunction with ^radar to obtain fine- 
scale wind measurements below 18 kn which is its n ^ximum altitude. The 
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purpose of the roughness elements is to control the random vortex shedding 
process which sends smooth balloons into erratic self-induced motion. The 
erratic aerodynamic, oscillations typically have periods of five seconds and 
amplitude: of a few ms * . The Jimsphere does not eliminate these oscilla- 
tions. Rather, they become more sinusoidal in character. As a result, it 
is easier to retrieve from wind profiles fine-scale structure of a scale 
equal to or less than the scale of the aerodynamical ly induced balloon 
oscillations. However, in the subcritical Reynolds number regime in the 
stratosphere the Jimsphere may offer little advantage over a smooth balloon. 

o Balloon Response 

Self-induced oscillations discussed above are the result of the interaction 
of a moving balloon with its own wake. It remains to consider the response 
of the balloon to variations in the environmental wind field through which 
it passes. Fichtl (1971) has made a linear perturbation analysis of the 
response of the Jimsphere to a vertically varying wind field. His analysis 
indicates that in the troposphere the Jimsphere is responsive to vertical 
variations in the wind field on the order of 10 m. Furthermore, Fichtl's 
analysis shows that the Jimsphere becomes less responsive at. higher altitudes. 
The analysis was not continued above 18 km and there appears to be no 
corresponding analysis for smooth balloons. Fichtl et al. (1972) have also 
considered the behavior of spherical balloons in wind shear layers and 
concluded that the Jimsphere may suffer a horizontal 11 velocity defect" as 
large as .6 m in the lower stratosphere. 

The response of tetroons to changes in the environmental wind field has 
been considered by several authors. The primary factor limiting the response 
of the tetroon to air motion is the restoring force of the balloon when it 
is carried away from its level of static equilibrium. Hanna and Hoecker 
(1971) have considered the response of constant-density balloons to sinusoidal 
variations in the horizontal plane of vertical wind speeds. Because of the 
restoring force, the balloon's vertical motion is less than the vertical 
motion of the air. More significantly, for sinusoidal motion the phase of 
the tetroon oscillation leads the phase of the air motion. These effects are 
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most pronounced for low frequency oscillations. For example, a sinusoidal 

-1 3 

vertical air motion with an amplitude of .5ms and a period of 10 sec 
will cause the balloon oscillation to lead the sinusoidal air motion by 30° 
and the vertical velocity of the tetroon to be only 57. less than the vertical 
velocity of the air. These results suggest that the analysis of flux data 
from turbulence statistics derived from tetroon data should be viewed with 
caution for periods greater than 10 or 15 minutes. Although horizontal 
velocities of the tetroon should be fairly representative of air motions, 
vertical velocities may be slightly underestimated and correlations between 
horizontal and vertical velocities may be poorly estimated. It is commonly 
assumed that horizontal and vertical velocity fluctuations measured from 
tetroon data are representative of air motions with periods of order several 
minutes or less. However, it should be noted that the sensitivity of the 
tetroon to the smallest scales of turbulence will decrease with altitude as 
density decreases. 

A. 3 DIFFUSION COEFFICIENTS FROM FLUX-GRADIENT RELATION 

If measurements of momentum flux are available simultaneously with a 
vertical wind profile then it is possible to estimate diffusion coefficients 
from flux gradient relations. According to this approach the vertical 
diffusion coefficient is defined by 



where u “ U-U and w = W-W are the departures from the mean of the longitudinal 
and vertical components of the wind. 

The vertical eddy momentum flux can be determined from tetroon data and the 

gradient of the mean wind can be evaluated from any wind profile data 

available. One possibility for deducing K from a single balloon flight 

z 

is to use the data from a tetroon approaching its "constant 11 flight level. 
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4.4 DIFFUSION COEFFICIENTS FROM LAG RANG IAN VELOCITY VARIANCE STATISTICS 

According to the theory of atmospheric diffusion as developed by 
Taylor (1921) the spread of smoke or some other passive additive can be 
related to the variance of turbulent atmospheric winds. Thus, in one 
dimension 

o 2 = U 2 ! 2 (4.2) 

for small diffusion time T where o is the half width of the cloud. If we 
define a one-dimensional diffusion coefficient K by the relation 

o 2 = 2KT , (4.3) 


then for small diffusion times 


— T 
2 A 


and it is possible to determine three-dimensional diffusion coefficients 
from the variance of measured velocity fluctuations 


*x - r 


(4.5a) 


Ky = v2 T 
1 2 


(4.5b) 


K, = w T 

Z T 


(4.5c) 


where u = U - U, v = V - V, and v e W - W are the departures from the mean 
of longitudinal, lateral and vertical velocities. A more general formulation 
of the one -dimensional diffusion is given by Taylor (1921) as 


2u2 a n v T > dTdt 


(4.6) 



where R x (t) is the Lagrangian auto-correlation function defined by 

R (T) + (4.7) 

2 

u . 

Clearly, the Lagrangian autocorrelation function will have its maximum 
value for t « 0 and will decrease with lag t as the velocities become 
uncorrelated. If R(t) ~ 1, as is the case for small I, Equation 4.6 
leads to 4.2, 

For general T it is useful to define 

I(T> - J* R/(T)dT (4.8) 

which has a value in the range 0 < I(t) < T. I(T) is dependent upon the 
structure of the Lagrangian turbulence spectrum. When the diffusion time 
exceeds the time scales for which there is any appreciable turbulent 
kinetic energy for a particular velocity component, the value of I(T) 
approaches an asymptotic value I* which is an integral diffusion time scale. 
In general this diffusion time scale will be different for each velocity 
component. 

If, following Taylor, the turbulence is assumed to represent a station- 


ary random process, then 



B . 

n 2 Ijj (T) 

(4.9a) 

Ky ", 

v 2 ly (T) 

(4.9b) 

and *Z 

w 2 ^ (T) 

(4.9c) 


define three-dimensional diffusion coefficients for arbitrary diffusion time 
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In the limit that the diffusion time becomes much longer than the time scale 
of the Lagrangian correlation function then the diffusion coefficients will 
approach asymptotic values 

= u 2 l£ (4.10a) 

- v 2 (4.10b) 

and Kg ■ w 2 I* for large T. (4,10c) 

“ ~ ~ 

Moreover, since w will be much less than either u or v , it is not surpris- 
ing that Kg values are many orders of magnitude less than or for 
global diffusion problems. 

4.5 DISCUSSION 

In this chapter a review has been presented of the application of 
balloon techniques to the measurement of small-scale turbulence and diffu- 
sion in the atmosphere. Several different approaches to tracking rising 
and floating balloons have been discussed and several problems involving 
the aerodynamics of rising and floating balloons have been considered. Two 
techniques have been presented for deri\ !.ng quantitative diffusion coeffi- 
cients from measured atmospheric winds. 

Until recently the most accurate technique for observing detailed 
atmospheric wind fields has been to track a Jimsphere with a radar such 
as the FPS-16. Endlich and Davies (1967) have explored the feasibility of 
using the FPS-16, Jimsphere system to measure turbulence in the free atmosphere. 
They concluded that the radar data showed sufficient resolution to observe 
fluctuating winds associated with light -moderate turbulence. No quantitative 
technique was presented to obtain turbulence parameters or diffusion 
coefficients from tracking ascending Jimspheres. 
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Any technique which employs vertically rising balloons suffers a 
major difficulty. Rising balloons simply do not stay in the same volume 
of air long enough to obtain representative turbulence statistics. To 
be sure, it is always possible to analyze the detailed structure of the 
vertical profile measured by the rising balloon. Nevertheless, it is 
extremely difficult to interpret the product of such an analysis. 

The most promising technique for obtaining quantitative turbulence 
measurements and diffusion coefficients involves the tracking of constant- 
level balloons. The floating constant -level balloon is carried in a trajec- 
tory which approximates the flow of the mean wind. High -resolution tracking 
data enables the separation of fluctuating longitudinal, transverse and 
vertical velocity components from mean values. Diffusion coefficients 
can be obtained from an analysis of Lagrangian turbulence statistics as 
outlined in Section 4.4, 

Velocities of less than 1 m must be observable. This precision is 
beyond the capability of most balloon-tracking systems discussed in this 
chapter without averaging over time periods of a minute or longer. Averag- 
ing over these periods is marginal for resolving turbulence of the vertical 
wind whose primary components include time scales of a few minutes. 




5. METRAC™ BALLOON TRACKING SYSTEM 

5.1 DESCRIPTION OF THE METRAC SYSTEM 

The METRAC system is based on the Doppler shift of a moving 
transmitter. Although the physical principles are well known, only the 
recent availability of low-cost digital components and UHF-VHF transistors 
have permitted an economically feasible electronic design. The METRAC 
system uses omnidirectional antennas for both transmitting and receiving 
and does not require mechanical or electronic scanning. This eliminates 
the elaborate pedestal and drive assemblies associated with dish antenna 
tracking systems. 

The basic elements of the METRAC system are an airtime trans- 
mitter and several receiving stations having known positions. As the trans- 
mitter moves in space, the frequency at each receiver equals the transmitted 
frequency plus a Doppler frequency shift, which is a linear function of the 
velocity of the transmitter. Because the true transmitted frequency may 
not be known, this Doppler shift cannot be determined from only the data 
at one receiver. However, the data from any pair of receivers permits a 
determination of the difference of received frequencies. Since these 
receivers are at rest with respect to each other, the frequency difference, 
equals the difference of the Doppler shift associated with the receiver ; 
pair. This Doppler difference is the only data required to determine the 
transmitter position relative to the receivers. 

The integrated Doppler difference associated with each pair of 
receivers is directly proportional to the slant range difference from the 
transmitter to each receiver, A known slant range difference determines a 
hyperbolic line of position on which the transmitter is located. The receivers 
are the foci of this hyperboloid. The data from three independent receiver 
pairs (four receivers) determines the transmitter position in space. 
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The electronics required to implement the system consists of an 
inexpensive balloon-borne transmitter, four or more receivers, and a central 
command console. The command console is used to record the Doppler data 
to determine the transmitter position. The METRAC system using six receivers 
is shown schematically in Figure 3. All receiver coordinates must be known 
accurately. An audio frequency communication link is used between the central 
command site and all receivers. A stationary radio transmitter generates 
a reference frequency; its coordinates are unimportant except that it must 
have a direct line of sight to each receiver. 

The balloon transmitter frequency, f , is nominally 403 MHz, which 
was chosen because this frequency falls in the band allocated for meteorolog- 
ical aids. The reference frequency, * re f> is kept about 2 KHz different 
from the balloon frequency by means of feedback from one of the receivers. 
Radiation from both transmitters is used by each receiver to form an output 
signal whose frequency f^^ is given by 


f roi 


£ t + 4 i - f ref 


(5.1) 


where d^ is the Doppler shift observed at this i t ^ 1 receiver. The signal 
to noise ratio is sufficiently large to permit multiplying this difference 
frequency by eight, thereby increasing resolution. The frequency is counted, 
repetitively sampled and stored on magnetic tape at the command console. 

The Doppler difference for any receiver pair yields 


Af - f - f . 
roi roj 


(5.2) 


The difference in counts given by 5.2 is used by the computer to solve for 
balloon position at the time of the samples, assuming that the intital 
location of the balloon is known accurately. 


HALLOON PACKAGE 



Figure 3. Schematic representation of the METRAC system deployment 










5.2 THE METRAC SYSTEM ACCURACY 

The Doppler differences obtained by forming the differences of 
frequencies received by pairs of receivers are the input to the METRAC 
solution algorithm. The transmitter location is obtained from this algorithm 
in a fashion similar to the solution of hyperbolic equations employed in 
common navigational problems. The potential accuracy of the METRAC system is 
determined by the effective wavelength of the transmitted signal and the 
geometry of the receiver array. 


/stem Geometry 


The basic measurement unit of the METRAC system is the wavelength, 

X, of the transmitted signal. Consequently, the maximum number of distinguish- 
able hyperbolic surfaces between two receivers separated by a distance D is 
D/X, and the transmitter position is determined to be on one of these surfaces. 
Increased resolution can be achieved by multiplying the received frequency 
by a factor M, The basic measurement unit is then X/M, denoted by X^, and 
the system accuracy increases as M increases. For a physically realizable 
multiplication factor of 8, X^ corresponds to a resolvable distance of 
approximately 10 cm at 403 MHz. 


As can be seen in the twordimensional example illustrated in 
Figure 4, the spacing between hyperbolic shells depends upon geometrical 
position. Only on the line between each pair of receivers is the separa- 
tion of the shells as small as X^, Clearly, the accuracy of the computed 
position depends upon the spacing of the shells as well as the orientation 
of the foci (receivers) with respect to each other and to the position of 
the transmitter. Maximum practical three-dimensional resolution is achieved 
by deploying an equally spaced ring of receivers about a central receiver. 

A minimum receiver array should consist of a triangle of receivers with a 
fourth in the center. The centrally located station is very important in 
giving good vertical resolution. When there are more than a minimum set 
of four receivers operating properly over some time interval, some of the 
solutions will be better than others because of the differences in the geometry. 






Figure 4. Two-dimensional examples of position uncer- 
tainty as determined by the intersection of 
hyperbolic shells. 
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The two-dimensional example of Figure 4 also shows that the hyper- 
bolic shells tend to become more nearly parallel to one another well outside 
the borders of the receiver array. This implies that the uncertainty in 
computing position of the transmitter grows as the balloon drifts away from 
the array. The uncertainty is largest in the radial direction. Increasing 
the averaging times in determining the winds when the balloon is well outside 
the station array will decrease the effect of the position uncertainty. 

e Sources of Error 

1 

There two primary sources of error which can affect the METRAC 
solution: count er.rors and station location errors. Count errors can come 

from a variety of sources, some interrelated. The most serious errors would 
arise if either the balloon or reference transmitter were interrupted even 
momentarily. This would cause dropped counts at all receiver stations 
simultaneously so that no solution could be computed. 

Count errors will also occur if the frequency difference (reference 
and balloon transmitter) tracking filter loses lock at any receiver. This 
will happen whenever the signal to noise ratio at the receiver becomes 
sufficiently small for an appreciable part of the time constant of the filter. 
Experience during field tests of a prototype METRAC system showed that this 
situation occurred most often when the transmitter was high and nearly 

t 

directly above a receiver station or when the transmitter was far (30-100 km) 
from the receiver. If the signal to noise ratio becomes small for a very 
short period of time, errors may occur in the counting even though the filter 
track remains locked to the frequency difference. Tests performed with a 
static transmitter showed that this type of error was generally random with 
l magnitude of only one or two counts. Errors in Doppler counts can also 

occur due to sampling uncertainty. However, these errors are non-accumulative 

•• and are also on the order of only one or two counts per sample. 
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An error from any of these sources will locate the balloon trans- 
mitter between the wrong hyperbolic shells. If the error is not random* as 
is the case for the first two errors described above, the position error 
will grow as the balloon gets further away from the array and the hyperbolic 
shells, get further apart. Random errors due to sampling uncertainty or 
occasional count errors will add artificial variance into the true position, 
but unless the balloon is well outside the baseline where the shells are 
far apart, the corresponding error in velocity is small even when computed 
over very short time intervals. 

In addition to position errors due to inaccuracies in the count 
of the Doppler frequency difference, errors also arise from uncertainties 
in the locations of the receiver stations themselves . The solution to the 
METRAC positioning problem requires frequency counts at a minimum of four 
receivers as well as initial launch coordinates relative to the receiver 
array. Because of the extreme resolution inherent to the tracking system, 
an error of ten meters in. the location of one station relative to the rest 
can be equivalent to as much as a 100 count sampling error* As the balloon 
moves away from its starting location, the position error will grow because 
the hyperbolic shells become more widely spaced as was discussed above. 

5.3 SUMMARY OF METRAC FIELD TEST RESULTS 

Field tests of the METRAC prototype system were carried out during 
the spring and summer of 1974. Figure 5 shows the locations on a map of 
the Twin Cities of the reference transmitter. (X) and the seven ground 
receivers (R) as they were deployed for these tests. The reference trans- 
mitter was installed on top of a 780 foot building in downtown Minneapolis. 
The receivers were located in a variety of commercial and residental 
neighborhoods throughout the Twin City area and were connected to the command 
site by leased telephone lines. 




represents the location of a receiver, X represents 
the location of the reference transmitter. 





Tests utilizing a stationary transmitter operating for extended 
periods of time demonstrated that improper counting of frequency cycles or 
"cycle slippage 11 occurred only intermittently as would be expected from 
the discussion of the previous section and was easily detectable. Utilizing 
more than the minimum number of four receivers makes it possible to correct 
for this type of intermittent problem. 

Figure 6 presents the results of an early test where a transmitter 
attached to a pole was carried along a prescribed path (dotted line) on the 
penthouse roof of a suburban hotel. Each one-second data point is plotted. 
The apparent systematic departure of a quarter meter may be due to the pole 
not being held vertically or more likely to an error in initialization. 

The "noise 11 in the track consists of both system errors and the wobble of 
the hand held pole. The si2e of these random errors agree very favorably 
with computer simulations of system errors. 

During March and April, 1974, nine test flights of the METRAC 
system transmitter were made. Eight of these consisted of radiosonde 
comparison flights in which both the METRAC system transmitter and a standard 
1680 MHz radiosonde were attached to the same balloon. The radiosondes were 
all tracked with a WeatherMeasure RD-65 tracking system borrowed from the 
University of Wisconsin. Many of these flights were also tracked with an 
optical theodolite. Ten additional test flights were made during June, 

1974. Results from several flights have recently been published (Gage and 
Jasperson, 1974) and only a sample will be discussed in this report. 

Figure 7 presents a comparison between 60 second METRAC and radio- ~ 
sonde winds (u, west-east component; v, south-north component) for the 
first 30 minutes of flight MF5. After an initial drift towards the ENE, 
the balloon traveled southward until, after 30 minutes, the balloon was 
approximately 10 km outside of the receiver array. The comparison between 
the radiosonde and the METRAC system derived winds are good below a height 
of 7 km after which the radiosonde winds tend to oscillate about the METRAC 




Figure 6. Trajectory obtained from KETRAC system data by moving. a METRAC 
transmitter along a prescribed path (dashed line) on the roof 
of a suburban hotel. 
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system winds. Flight MF5 was also tracked with a theodolite and Figure 8 
shows the comparisons for 20 second winds for a section of the.f light. 

This section was chosen because of the particularly large variation found 
in the theodolite winds. However, the METRAC system winds confirmed the 
erratic wind structure except for one bad theodolite reading at 980 seconds. 
Because both systems were tracking the same balloon, this analysis, of 
course, says nothing about the time or horizontal space scale over which 
this large amount of wind speed structure existed. Figure 9 presents the 
wind profile for the same flight up to 20 kilometers in height. No 
comparison data is shown in this figure because of the condensed scale 
necessary to illustrate the detail present throughout the entire profile. 

A considerable amount of sharp wind velocity structure clearly exists above 
the tropopause height of 10,8 kilometers ; 

The high resolution capability of the METRAC system is further 
illustrated in Figure 10 which shows each onersecond measure of balloon 
borne transmitter velocity for a one -minute dat* sample. The regular 
oscillation seen in this figure shows the pendulum motion of the transmitter 
which was suspended several meters below the balloon. 













6. DIFFUSION COEFFICIENTS FROM A SAMPLE OF METRAC TETROON FLIGHT DATA 

During June 1974 two constant-level balloons were launched in Minneap- 
olis and tracked using the Minneapolis METRAC system. In this section the 
results of an analysis of Lagrangian turbulence statistics from a ten 
minute segment of MF18 are presented. The data used in this analysis are 
reproduced in Figure 11. At the top of Figure 11 are plotted 10 second 
samples of the vertical velocity between 700 and 1300 seconds after launch 
of the tetroon. The mean altitude of the tetroon during this portion of 
flight was about 1.75 km above the surface. In the middle of the figure 
are plotted 10-second samples of the lateral component of the fluctuating 
velocity. At the bottom of the figure the longitudinal component of the 
fluctuating velocity is reproduced. 

6.1 DIFFUSION COEFFICIENT FROM FLUX -GRADIENT METHOD 

The turbulent velocity data shown in Figure 11 can be used to 

obtain momentum fluxes. The result is plotted in Figure 12. The vertical 

4 

momentum flux uw for the period of interest is approximately ,03 x 10 

2-2 4 

cm s . During this segment of flight the tetroon rose 3 x 10 cm and 

2 -1 

the longitudinal component of velocity decreased roughly 10 cm s 
Therefore, an estimate of from the flux gradient method is given by 

K« ~ .9 x 10 5 cm 2 s* 1 . 

M 

AZ 

There is no way to estimate and Ky from the available flux data since 
horizontal gradients were not measured. 

6.2 DIFFUSION COEFFICIENTS FROM ANALYSIS OF VARIANCE STATISTICS 

The Lagrangian velocity can be analyzed in several ways. A 
common approach is to form the Lagrangian autocorrelation function defined 
by Equation 4.7. This has been done for the first five minutes of the 
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Figure 11. Departures from the mean of vertical, 
longitudinal and transverse velocities 
as measured from 600 seconds of tetroon 
flight MF18 launched at 1605 CDT on 
13 June 1974. 
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Figure 12. Horizontal and vertical momentum fluxes 
calculated from data presented in Figure 
11 . 
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segment of tetroon flight MF18 reproduced in Figure 11 and the results are 
shown in Figure 13. In this figure the normalized correlation function is 
plotted with lags ranging from 0 to 250 seconds. The correlation coefficient 
for the longitudinal velocity shows evidence of a periodic structure with a 
period of about 160 seconds. The transverse and vertical velocity correla- 
tion functions appear qualitatively similar with negative correlations for 

lag greater than 130 seconds. The variance of the 10-sec samples of fluctu- 

4 2 2 

ating velocities were .13 x HT ctn s~ for the vertical component, .63 x 

4 2 —2 4 2 — * 

10^ cm s~ for the transverse component and .58 x l(j cm s for the 

longitudinal component. These values can be used together with estimates 

of I(T) to obtain estimates of the diffusion coefficients according to 

Equation 4.9. The results are summarized in Table 1. 


TABLE 1. 

DIRECTION VARIANCE 

Longitudinal .58 x lO^cm 2 s~ 2 

(X) 

Lateral .63 x lO^cm 2 s~ 2 

00 

4 2-2 

Vertical .13 x 10 cm s 

(Z> 


I(T) K 

30 sec 1.8 x lO^cm 2 s~* 


50 sec 3.2 x lO^cm 2 s~* 


60 sec 


.78 x 10 5 cm 2 


-1 

s 
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6.3 DISCUSSION 

The high resolution three-dimensional positioning capability of 
the METRAC tracking system permits meaningful velocity variance data to be 
collected on time scales as small as a few seconds. Figure 14 presents a 
computer derived maximum random positioning error distribution associated 
with system frequency count, errors of magnitude 1. Maximum positioning 
errors in each of the three components x (E-W), y (N-S) and z are shown 
with respect to the four Minneapolis field test stations (dots with lines 
connecting the perimeter) used to compute the tetroon statistics. The 
error magnitudes are valid for a one kilometer height. The arrow in the 
bottom third of each sub-figure shows the horizontal trajectory of the 
tetroon between 700-1300 seconds into the flight. 

It is clear from Figure 14 that random positioning errors in x and 
y associated with single count errors are negligibly small. In reality as 
was discussed in Chapter 5, count errors are often one count or less but 
can through both round-off and sampling errors reach two counts. This 
implies maximum random errors of two meters. in the horizontal and about 
six meters in the vertical. For 10-second sampling interval, this yields 
error bounds of 0.2 ms^in the horizontal and 0.6 m s* 1 in the vertical. 

The average or expected value of this error has been determined to be about 
a factor of five less than the maximum errors, thus yielding expected 
velocity errors in all three components of less than 0.1 ms* 1 for the 
10-second sampling intervals used in the preceding analysis. The variances 
in Table 1 should be in error less than 17*.. The uncertainty in the diffusion 
coefficients given there should be dominated by the uncertainty of determin- 
ing 

Si R(T)dr. ' . ' 

The agreement between the vertical diffusion coefficient obtained from the 
flux-gradient method and the vertical diffusion coefficient obtained from 

• i . .. • . ' ; ■ v . 
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the Lagrangian turbulence Is probably coincidental. Much more confidence 
can be attached to the estimate based on the Lagrangian turbulence statis- 
tics since it is undoubtedly a better measure of representative atmospheric 
turbulent structure. Nevertheless , a much greater sample of data would 
need to be analyzed in order to provide a statistically representative value 
of the three-dimensional diffusion coefficients which are typical for the 
altitude, latitude and season for which this sample was recorded. 
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7. APPLICATION OF THE METRAC SYSTEM FOR TRACKING CONSTANT- LEVEL BALLOON S 
IN THE STRATOSPHERE 


The purpose of this chapter is to discuss the considerations pertinent 
to utilizing the METRAC balloon-tracking system to collect stratospheric 
turbulence data from tetroon flights, A primary consideration is to assure 
a high probability of a successful flight. Also, accuracy is desired at 
stratospheric heights. And finally, it is important to strive to achieve 
the maximum length of tracking time for each flight. 

As was discussed in Chapter 5, a minimum number of four METRAC system 
receivers must operate properly over each sampling interval in order to 
compute a new position solution. If there are less than four receivers 
operating over that interval, the new position must be extrapolated from 
previous positions. If the sampling interval is any appreciable length 
of time, say even 10 seconds, extrapolation can be dangerous because all 
future position calculations depend upon the previous computed or assumed 
positions. During the Minneapolis field tests the problem of momentary 
receiver losses was eliminated by using seven receivers to track the 
balloons. Utilizing more than the minimum number of four receivers would 
also be necessary to assure a high probability of success for stratospheric 
tetroon flights. 

The METRAC positioning system has an advantage over all of the tradi- 
tional angle tracking systems in that velocity accuracy is not strictly a 
function of distance from the antenna and elevation angle. With the METRAC 
system the resolution is dependent upon the position with respect to the 
entire array of receivers. The resulting errors or error volumes may be 
described in many ways. Figures 15-18 present the positioning errors as 
a function of geometrical location with respect to a triangular array of 
receivers represented by open circles. They are expected or average errors 
associated with each full Doppler count of error. Only the errors in the 
x (E-W) and z directions are illustrated. Since the station array is 
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Figure 16. Horizontal distribution. of expected random 
positioning error per Doppler count error 
in x (E-*W) and z at a height of 10 km. Th 
open circles represent the stationary 
receiver array. 
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Figure 17. Horizontal distribution of expected random 

positioning error per Doppler count error in 
x (E-W) and z at a height of 20 km. The open 
circles represent the stationary receiver 
arrav . 
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Figure 18. Horizontal distribution of expected random 

positioning error per Doppler count error in 
x (E-W) and z at a height of 30 km. The open 
circles represent the stationary receiver 
arrav . 
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•symmetric, the y (N-S) errors are similar to the x errors. 

Figures 15-18 show that the horizontal errors in positioning with 
four receivers configured in a triangular array remain reasonably consis- 
tent with height and are small at large distances from the array itself. 

The height errors are interesting in that they are largest at low levels 
outside of the receiver array. At heights of 20 to 30 km these positioning 
errors are less than two meters at distances several tens of kilometers 
away from the receiver array. 

Several factors must be kept in mind in order to discuss maximizing 
the length of useful data for each flight. These factors will be discussed 
in the following paragraphs. 

When discussing stratospheric turbulence, the altitude range under 
consideration can extend from 7. to 50 km. The lower limit in particular 
will depend upon the season and latitude. However, even for the lower 
stratospheric heights, considerable time will be required to let a normal 
super-pressure tetroon rise to the desired altitude. Two approaches have 
been used in launching tetroons. These include towing the constant-level 
balloon to the desirable height with a more buoyant balloon and/or releasing 
the tetroon upwind so that it is at or near float altitude by the time it 
approaches the radar site. Either of these two techniques can also be: used 
in conjunction with the METRAC system. In the case of the upwind launch, 
the original launch coordinates must be known and the launch point must 
be within radio line-of-sight to at least four receivers. This last 
requirement can be more restrictive than one might desire. Some primitive 
computer simulations have shown that it may be possible to circumvent know- 
ledge of the original launch location and the launch line-of-sight require- 
ments altogether if more than four receivers operate properly for a period 
of time after launch. However, this has not been verified with field 
.experiment data. 

Another ^factor which must be considered in gathering stratospheric 
turbulence data with the METRAC system is that the transmitter must be 
able to operate without substantial power loss for extended periods of 
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time (1-3 hours) at very cold temperatures (-40 to -70°C). Special care 
should be maintained in temperature-compensating the transmitter circuitry, 
thermally insulating the transmitter from ambient air temperature and. 
providing a chemical heat source for the battery power supply. 

A final but very important factor which must be considered in using 
the METRAC system to collect stratospheric data is the receiver array 
design. Actually the design is determined by requirements already stated. 

In summary these include: 

1) Redundant (more than 4) receivers 

2) Triangular arrays with interior stations 

3) Accurate data over long trajectories 

In practice there is an implied fourth requirement: To minimize the cost 

of operating, maintaining and processing data from the system (i.e. minimize 
the number of receivers). Figure 19 presents a surface receiver station 
configuration which fulfills all four of these requirements. A sacrifice 
which is made in this configuration is that it is not symmetric. For best 
results the array should be oriented along the direction of the mean winds 
at the levels of interest. Accurate wind measurements over distances in 
excess of 100 km should be possible. This will yield useful statistics in 
Lagrangian turbulence on time scales as short as ten seconds over flight 
periods of one-half to three hours, depending upon the wind speed. 
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8. CONCLUSIONS AND RECOMMENDATIONS 

This report has presented a review of the basis for present estimates 
of diffusion coefficients at stratospheric altitudes. Horizontal diffusion 
coefficients are dominated by large-scale eddy exchange processes and can 
be evaluated from available meteorological data. However, vertical diffusion 
coefficients are dominated by small-scale turbulent diffusion. The turbulence 
responsible for vertical diffusion coefficients has not been adequately 
sampled at stratospheric altitudes and only crude estimates exist of global 
mean vertical diffusion coefficients based on the behavior of radioactive 
tracers. A new technique to measure small-scale turbulence and diffusion 
at stratospheric altitudes is required. 

From a review of alternative techniques for obtaining turbulence 
measurements it is concluded that high-resolution tracking of constant- 
level balloons at stratospheric altitudes could provide a simple, direct way 
of measuring Lagrangian turbulence statistics. The Lagrangian velocity 
variance statistics lead directly to a quantitative measure of three- 
dimensional diffusion coefficients. 

The high-resolution capability of the METRAC system for tracking 
tetroons in the troposphere has already been demonstrated and a sample of 
the tetroon data has been included in this report. An error analysis of 
the positioning capability of the METRAC balloon-tracking system leads to 
the conclusion that this system will yield sufficient accuracy to measure 
turbulence at stratospheric altitudes even for data averaged over as little 
as ten seconds. 

It is concluded, therefore, that the METRAC system can be used to obtain 
three-dimensional small-scale diffusion coefficients at stratospheric alti- 
tudes. It is recommended that a field demonstration be undertake.! to use 
METRAC to collect samples of stratospheric turbulence data in order tc 



compute three-dimensional diffusion coefficients* If the results of the 
field demonstration are successful, the METRAC system should be employed 
to directly measure stratospheric turbulence at selected geographical 
locations and seasons in order to build up a climatology of the variability 
of diffusion coefficients* 
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